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THE  PROBABILITY  OF  MULTIPLE  CORRECT  PACKET  RECEPTIONS 
IN  A  ML^LTIRECEIVER  FREQUENCY-HOPPED 
SPREAD-SPECTRUM  SYSTEM 

1.  INTRODUCTION  AND  PROBLEM  FORMULATION 

An  important  quantity  in  spread-spectrum  radio  networks  is  the  probability  that 
exactly  /  out  of  m  packet  transmissions  8ire  successful,  given  that  k  users  attempt  to 
transmit  their  packets  simultaneously;  this  quantity  is  denoted  by  P(*’,m  —  l\k).  This 
quantity  is  essential  for  the  integration  of  the  first  three  layers  (physical,  data  link,  and 
network  layer)  of  the  ISO  network  model  for  packet  radio  networks  using  spread-spectrum 
signaling  and  forward-error-control  and,  as  such,  it  enables  the  design  and  performance 
evaluation  of  multiple-access  protocols  for  such  networks. 

The  integer  m  in  P{l,m  —  /|fc)  denotes  the  number  of  receivers  of  interest;  in  most 
practical  situations,  m  <  k.  Specifically,  in  problems  involving  multireception  with  a  bank 
of  m  receivers  at  a  single  location,  the  probability  mass  function  (pmf ) 

P{l^m  —  l\k)  for  /  — 0,  and  m<k 

describes  the  multireceiver  performance.  Moreover,  in  problems  in  which  the  evaluation 
of  the  throughput  or  delay  of  various  packet  radio  network  protocols  is  desirable, 

is  required,  where  Pc{l\k)  denotes  the  probability  of  tmy  I  correct  packet  receptions  out  of 
k  simultemeous  transmissions  (see  [1],[2]). 

Consequently,  in  practical  spread-spectrum  packet  radio  networks,  there  is  an  undis¬ 
puted  need  to  evaluef#>  the  probabilities  P(/,  m  — /|A:)  and  P7’(/,  m  — /U  )  =  ^  ^  ^P(/,  m  — 

l\k)  (/  =  0, 1,  ■  •  •  ,m  cmd  m  <  k)  for  different  spreading  signaling  formats,  data  modulation 
schemes,  and  error-control  coding  schemes. 

In  this  report,  we  evaluate  these  quantities  for  frequency-hopped  (FH)  spread-spectrum 
multiple-access  (SSMA)  networks.  Specifically,  we  develop  an  exact  expression  for  P(/,  m  - 
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^P{l,k-l\k)  for  l  =  0,l.-  -.k 
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l\k)  eind  aii  approximation  based  or  Gaussian  multivariate  densities.  The  exact  expres¬ 
sion  is  difficult  to  compute  and  its  computational  complexity  grows  exponentially  with 
m.  By  contrast,  the  Gaussian  approximation  is  computationally  efficient  and  its  com¬ 
plexity  grows  linearly  with  m.  Numerical  results  obtained  from  these  two  methods  are 
compared  with  those  obtained  via  the  independence  assumption  method  commonly 
used  in  the  literature  [1],[2].  This  method  assumes  that  packet  errors  among  different 
receivers  are  mutually  independent,  which  greatly  simplifies  the  computation.  We  further 
establish  that  the  independence  assumption  can  be  trusted  in  most  cases  and  derive  the 
range  of  parameters  in  which  each  of  the  two  approximations  (Gaussian  and  independence) 
is  preferable. 

Derivations  and  comparisons  are  carried  out  for  FH/SS  systems  employing  MFSK 
modulation  with  noncoherent  demodulation  and  Reed- Solomon  (RS)  {n.kc)  forward  error- 
control  coding  with  erasures-only,  errors-only,  and  errors  /erasures  minimum-distance  [3] 
decoding.  It  is  assumed  that  each  RS  symbol  carries  one  M-axy  symbol  (i.e.,  n  =  M),  that 
each  FH  dwell  time  (hop)  carries  one  RS  symbol,  and  that  one  RS  codeword  per  packet  is 
transmitted.  The  frequency- hopping  patterns  of  the  different  users  are  modeled  as  random 
memoryless  hopping  patterns  [4].  Thus,  each  of  q  available  frequencies  are  visited  with 
equal  probability  emd  independently  of  each  other  during  any  dwell  time  (hop)  by  each  user 
and  mutually  independent  hopping  patterns  are  assigned  to  distinct  users.  The  \'arious 
users  are  packet -synchronous  but  may  be  hop-asynchronous;  in  this  context,  both  hop- 
synchronous  cuid  hop- asynchronous  FH/SSMA  systems  are  considered.  Also  thermal  noise 
modeled  as  additive  white  Gaussian  noise  (AWGN)  is  incorporated  in  the  analysis. 

This  report  is  organized  as  follows;  In  Section  2  exact  expressions  for  P(/,m  —  l\k) 
are  derived  for  all  cases  of  interest  enumerated  above.  In  Section  3,  the  corresponding 
expressions  beised  on  the  Gaussian  approximation  technique  are  derived.  In  Section  4. 
the  approximation  based  on  the  independence  assumption  is  cited.  Num.erical  results  and 
comparisons  of  the  three  approaches  are  presented  provided  in  Section  5.  In  Section  C. 
conclusions  are  drawn. 
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2.  DERIVATION  OF  EXACT  MULTIRECEPTION  PROBABILITIES 


We  axe  interested  in  finding  an  exact  expression  for  the  probability  P{l,m  —  l\k)  of 
I  receivers  receiving  correctly  and  m  —  I  ones  receiving  erroneously,  for  a  specific  set  of 
receivers.  It  eissumed  that  aJl  users  employ  similar  MFSK  modulation  (with  M  frequency 
tones)  with  noncoherent  demodulation  and  identical  extended  RS  (n,  fcc)  codes  with  c<^de- 
word  length  n  symbols  {n  =  M)  and  kc  information  symbols  per  codeword.  (Refer  to  the 
additional  assumptions  in  the  previous  section.)  Due  to  the  symmetry  in  the  system,  we 
can  equiv2dently  find  the  probability  of  the  first  I  receivers  decoding  correctly,  while  the 
remaining  m  —  I  receivers  decode  in  error.  In  the  sequel,  we  implicitly  assume  that  m  >2. 
For  m  =  1,  the  model  reduces  to  the  single  receiver  model  analyzed  in  [4]. 

For  FH/SSMA  communications,  the  probability  of  a  coded  symbol  error  is  upper- 
bounded  by  the  probability  of  a  hit,  which  is  a  function  of  the  available  frequency  slots  q 
and  the  number  of  contending  users  fc,  where  k  >  \.  Subsequently,  we  denote  by  Ph{q.,k) 
the  probability  of  a  hit.  Recall  from  [4]  that  for  hop- synchronous  FH/SSMA  systems 

PUq,k}  =  l-il-l/g)^-^  (2.1a) 

and  for  hop- asynchronous  systems 


It  is  assumed  that  k  >  m  >2  and  q  >  m  >  2  so  that  the  above  quantities  are  guaranteed  to 
be  nonnegative.  For  8isynchronous  systems  the  above  upper  bound  is  valid  for  any  number 
of  symbols  per  hop  (dwell-time)  equal  or  larger  than  1  (slow  frequency-hopping). 

2.1.  Errors-Only  or  Erasures-Only  RS  Decoding 

The  ith  receiver  receives  correctly  the  transmitted  packet,  if  the  number  of  hits  h(i), 
for  1  <  i  <  rn,  satisfies 

0  <  h{i)  <  t, 


where  t  denotes  the  correction  capability  of  the  (n,kc)  block  code  {kc  is  the  number  of 
information  symbols  and  n  the  total  number  of  symbols  per  codeword)  For  pure  error- 
correctioii, 


t  = 


(2.2a) 
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while,  for  erasure  correction, 


t  =  djnin  -  I  =  n  -  kc.  {2.2b) 

In  the  synchronous  case,  each  receiver  output  depends  only  on  other  receivers  out¬ 
puts  during  the  same  dwell  time.  This  applies  also  to  the  asynchronous  case,  if  proper 
interleaving  takes  place.  Then  the  total  system  operation  becomes  memorjdess. 

As  amplified  in  Appendix  A, 

P{l,m-l\k)  = 


—  h{l)  <  t, . . .  ,0  <  h{l)  <  t.,t  +  1  <  h{l  1)  <  n, . . .  ,i  +  \  <  him)  <  n) 


=  EE  ■■  E  E 


n  \  /n-(i 

fJV 


2”'  -2 

n  -  ^nn 

nn  =  1 
^2'"-l 


■  {Po)  ""=» 


^  n 


where  P„n  =  P{Enn)  denotes  the  probability  of  the  event  Pnm  under  which  the  m  de¬ 
modulator  outputs  during  a  particular  symbol  of  the  codeword  (packet)  correspond  to  the 
binary  representation  of  nn  (recall  0  and  1  denote  correct  and  incorrect  reception,  respec¬ 
tively),  and  is  the  number  of  times  the  event  E„n  occurs  in  one  codeword.  Of  course, 
any  other  correspondence  of  the  above  events  and  the  natural  numbers  would  work  as  well. 
Note  that  in  (2.3),  all  events  having  same  weight  have  equal  probabilities,  although  this 
does  not  simplify  the  expression. 

The  range  of  for  the  sums  in  (2.3)  is  to  be  obtained  from  a  Diofantine  analysis  of 
the  inequalities; 


0  <  ^  ,  i  =  1,2,.../ 


(2.4a) 


t  -f-  1  <  ^  •  £„„  ,  i  =  /  -t-  1,  /  +  2, . . . 


i2.4b} 
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where  Qnn  =  1  or  0,  according  to  whether  the  ith  component  of  the  vector  event  E„n  is 
1  or  0,  that  is,  it  takes  part  in  the  ith  receiver  error  count  or  not.  We  have  to  add  another 
constraint  to  the  m  constraints  posed  by  (2.4a,2.4b),  namely  that 


0  <  ^nn  <  n 


(2.4c) 


The  purpose  of  this  constraint  is  to  ensure  that  we  do  not  surpass  the  codeword  length  n 
by  permitting  highei  values  of  the  ^nnS. 

It  remains  to  find  expressions  for  all  P{Enn)-,  for  nn  =  1, 2,  •  •  • ,  2'”.  This  is  equivalent 
to  finding  the  probability  of  having  p  demodulator  outputs  correct  and  m  —  p  ones  in 
error  during  a  particular  transmitted  symbol  (of  identical  order  for  all  receivers),  for  p  = 
0, 1,  •  •  •  ,  m.  These  probabilities  should  be  a  function  of  p,  m,  k.  and  q.  We  denote  them 
Ps[py'^,<},k).  First  we  find  P{Eq)  =  Pcc...c  —  0, g,  fc),  that  is.  the  probability  of 

decoding  correctly  all  the  simultaneous  symbols  in  all  receivers  ;  recall  Eq  corresponds  to 
the  vector  event  (0, 0,  •  •  • ,  0).  Because  of  the  symmetry  we  get 

Pc^^c  =  P{c\c^^c)  ■  Pc^^c 

m  Tn  —  1  m  —  I 

=  P(clc^^^^c)  •  P{c\c^^^c) . . .  P{c\cc)  ■  P(clc)  •  P(c) 

m  — 1  m— 2 

m 

=  +  +  (2.5) 


In  the  above  and  in  subsequent  expressions  the  following  notation  is  used;  P(c|c^^^^c)  or 

t 

P(e|c^^^^^c),  i  =  l,2,..,m  —  1,  denote  the  conditional  probabilities  of  a  particular  symbol 

1 

of  a  single  receiver  being  correct  or  incorrect  given  that  symbols  of  the  same  order  of  i 
other  receivers  axe  correct;  similarly  P(e,^^^e|c  . . .  c)  denotes  the  conditional  probability 

j  « 

of  a  particular  symbol  of  j  receivers  being  incorrect  given  that  symbols  of  the  same  order 
of  i  other  receivers  are  correct;  and  P(c  . . .  c)  or  P(e  . . .  e)  denote  the  the  unconditional 


5 


(absolute)  probabilities  of  a  particular  symbol  of  i  receivers  being  correct  or  incorrect, 
respectively.  We  now  find  P{E2’^-i)  —  Pee  -e  =  P-(0,m,q,k}.  We  have 

Pe...e  =  P(ele  .  e)P,  . .  e 


=  1-P(cie.  ..e)  Pe... 


,  P(e...e|c)-P(c)  „ 

_  i 

•*  ...  e 


-  P(...e  -  P{c]  ■  P{e. .  .t\c, 


=>  P3{0,m,q,i)  =  P^(Q,m  -  \,q,k)  -  [I  -  Ph{q,k)]  ■  Ps(Q,m  -  l,q  -  l,k  ~  1)  (2.6) 

In  these  expressions  Equation  (2.6)  is  a  recursive  formula  for  finding  P,(0,  m,  q,  fc).  The 
solution  of  this  equation,  as  shown  in  Appendix  A,  is 


m  /  \  ’ 

p,{o,m,q,k)  ^  I  +  Y.  +  + 

1=1  ^  ^  j=i 

Proceeding  one  step  further  we  obtain  the  more  general  expression 
P{Enn)  =  Pa(p,m  -  p,q,t) 


=  p 


C  .  .  ,  C  C  .  .  .  C 


=  P(e  . . .  ejc  . . .  c)Pc  . . . 


m-p  p 


=  Ps{0,m  -  p,q  -  p,k  -  p)  ■  Ps(p,0,q,k] 


=  Yl[l-Pk{q-J  +  l,k-j  +  l)] 


m-p  ,  ^  I 

•  1  +  E  ['^~^)Y[[l-Ph{q-p-3  +  l.k-p-J+  1)1 

.=  1  ^  ^  j=\ 

L  ”* 


6 


(2.8) 


The  equality  P(e^^^ejc^^^c)  =  Pj(0,  m  —  p,q  —  p,k  —  p)  is  obtciined  by  observing  that 

m-p  p 

conditioning  on  a  particular  symbol  of  p  receivers  being  correct  is  equivalent  to  reducing 
the  number  of  available  frequencies  for  hopping  q  and  the  number  of  users  k  by  p.  The  last 
equality  in  (2.8)  is  then  obtained  through  substitution  from  (2.5)  and  (2.7).  Equations  (2.3) 
and  (2.8),  together  with  the  constraints  posed  by  (2.4a),  (2.4b),  (2.4c),  give  the  solution  to 
our  problem.  As  it  becomes  cleax  from  (2.3),  the  exact  evaluation  of  P[l,m  —  l\k)  requires 
the  computation  of  2"*  —  1  dependent  sums,  in  which  the  limits  should  be  found  through  a 
Diofantine  analysis  of  (2. 4a), (2. 4b),  and  (2.4c).  In  addition,  the  summands  are  powers  of 
P(Pnn),  which  can  be  computed  through  (2.8).  Due  to  those  computational  requirements, 
exact  expressions  are  nearly  impossible  to  evaluate,  for  m  >  4.  However,  as  the  next 
sections  indicate,  useful  approximations  have  been  developed  with  satisfactory  accuracy  in 
diiferent  ranges  of  the  various  system  parameters  of  interest. 

2.2.  Inclusion  of  AWGN  and  Errors-Only  or  Erasures-Only  RS  Decoding 

The  expression  in  (2.3)  remains  valid  in  its  general  form;  however,  the  expressions  for 
Pe  ...  e  c  ...  c  change  in  order  to  take  the  effects  of  AWGN  into  account.  Let  us  find  first 


Pc  ...  c-  For  that  we  get 


Pc...c  =  (^-p/v)‘'-Pc'';l. 


(2.9) 


where  PcV..  c  denotes  the  previously  evaluated  Pc  ...  c  in  (2.5)  by  counting  errors  (or 

t' 

erasures)  from  hits  only,  and  P^  is  the  symbol  error  probability  of  the  MFSK  due  to 
AWGN.  To  find  P,  . .  ,  we  write 


Pt^^c  =  P(e|e^^e)  •  Pe^e 

^  M-i  ,.-1 

=  Pc  ...  c-  P{c)  ■  P(c^^e|c) 
^-1 


=>  P,{fi,q,k)  =  Pe{^x-l,q,k)-{l-  Pr,)[l-  PH{q,k)]-  Pe{n-l,q-l^~  -1).  (2.10) 
The  recursive  equation  (2.10)  is  similar  to  (2.6)  and  thus  has  as  solution  the  expression 


which  gives  the  desired  expression.  For  our  problem,  n  +  u  —  ni. 

The  probability  P\  of  MFSK  symbol  errors  due  to  AWGN  given  by 


(2.12) 


P,.  = 

^  \  m  J  Tn  + I 

fn:=  I  ' 


(2.13) 


where  Eb/No  is  the  information  symbol  signal-to-noise  ratio  and  r  is  the  code  rate  of  the 
system  (r  =  ^'c/n). 

2.3.  Errors  and  Erasures  RS  Decoding 

In  this  section,  we  analyze  the  case  in  which  all  receivers  employ  combined  errors  and 
erasures  decoding.  Thus  during  the  minimum  distance  RS  decoding  correction  of  both 
errors  and  erasures  is  attempted.  We  assume  *hat  erasures  happen  whenever  a  hit  from 
other  users  occurs  and  they  are  detected:  errors  are  caused  by  AWGN  only. 


S 


It  is  well-known  that  any  RS  code  can  decode  correctly  any  received  word  having  5 
erasures,  t  errors  if  these  numbers  satisfy 

2t  +  s<  -  1  =  n-k,  (2.14) 

We  start  the  analysis  for  this  case  by  noting  that  the  basic  “simultaneous’'  events  in  the 
multi-receiver  are  now  defined  by  x  =  (xi,...,Xm).  ''iiere  x,,  ?  =  can  take 

vzdues  c,  e,  s  corresponding  to  correct,  error,  or  erasure  symbol,  respectively.  Again, 
because  of  symmetry,  the  probability  of  finding  in  a  particular  symbol  (or  dwell-time) 
some  receivers  in  error,  some  in  correct  receptioa,  and  some  in  erasure  is  independent  of 
the  particular  order  of  the  receivers. 

For  a  given  symbol  we  define 

Pr(  Uf  receivers  in  error,  Ug  receivers  in  erasure,  Uc  receivers  correct  ) 


where,  of  course. 


n<  r»c 


P{ne,ng,nc) 


ne  +  Og  He  =  rn. 


(2.15) 


(2.16) 


To  find  P{!,m  —  l\k)  we  notice  that  the  basic  difference  from  the  treatment  in  Section  2.1, 
is  that  now  we  have  three  different  states  in  every  symbol.  So,  instead  of  having  2”*  —  1 
sums,  we  have  3"*  —  1  sums.  If,  we  define  as  E„„  the  event  according  to  which  the  m 
demodulator  outputs  during  the  saxne  symbol  correspond  to  the  ternary  representation  of 
nn,  where  0  denotes  correct  reception,  1  denotes  an  erasure,  and  2  denotes  an  error,  and 
set  to  denote  the  number  of  times  the  event  E„n  occurs  in  one  codeword,  then  we  can 
write 


P{1,  rn 


'w  =  EE-  E  E 

(\  ti  ^3"*  _  3  tarn  _  , 


n  — 
(2 


3^  -  1 


n  — 


(Po) 


E 
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(2.17) 


The  range  of  values  of  ^2,  •  ■  can  be  found  from  the  m  inequalities: 


0  <  ^  •  fnn  <  «  -  ^'c  for  2  =  1, 2,  •  ■  • ,  / 


(2.1Sa) 


n  -  4-  1  <  ^  ■  („n  for  i  =  /  +  1,  /  +  2,  •  ■  •  , : 


(2.186) 


where  anl  takes  values  G  {0,1,2},  depending  on  what  type  of  event  (correct  .symbol, 
erasure,  or  error)  is  implied  for  receiver  t  £  mj  from  the  z-th  component  of  the 

ternary  representation  of  7172  (identifying  the  event  £’nn);  that  is,  ul/n  =  0  for  correct 
reception,  =  1  for  erasures,  and  ai/J  =  2  for  errors.  To  the  above  inequalities  we  must 


add  the  condition 


0  <  ^  fnn  <  n 


(2.18c) 


which  ensures  that  the  total  number  of  symbols  remains  smaller  than  n,  the  codeword  (or 
packet)  length. 

In  order  to  evaluate  P{l,7y7  —  /)A’).  we  need  to  calculate  expressions  of  the  form 


e  ...  e  s  ...  s 


r  ...  Cl  where  71^  +  ri^  +  71^  —  m.  We  can  write 


n  g  n .  Of 


Pe  ...  fa  ...  a  c  ...  c  =  ■■  ■  s\e  .  .  ■€  C  .  .  c)  •  F(e  .  .  cjc  .  .  .  c)  •  P(c  .  .  .  C )  (2.19) 


Tig  n,  rif 


n.  rig  Hr 


rig  ric 


where 


P{C  C)  =  ^{(1  -  P^^.)[l  -  P,(g  -  j  +  IJ-  -J  ^  1)]} 


(2.20a] 


P(€  ...  e\c  c)  =  l[{PN[l-Ph(q-n,-j  +  l.k-72,-j+l)]}  (2.206) 


n.  n. 
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P{s...s\e..€  C...C) 


n ,  Tic 


=  1  +  ^(-1)' f”*")  j[[l  -  -  ne  -  nc  -  j  +  l,fc  -  ne  -  He  -i  +  1)] 


.=  1 


J=1 


(2.20c) 


so  that 


P{E^n) 


P(ne,n,,nc) 


j=i 


•  -  Uc  -  J  +  l,k -He  -  j  +  1)]} 

J=t 

■  1  +  ~  ■  j  +  l,k  -  Tie  -  Tic  -  j  +  1)] 


(2.21) 

2.4.  Asynchronous  FH/SSMA  Case 

It  is  straightforwcLrd  to  extend  the  results  of  the  hop-synchronous  ceise  to  the  hop- 
asynchronous  one,  if  we  assume  that  all  symbols  wnthin  each  codeword  are  interleaved.  In 
the  asynchronous  case,  in  addition  to  the  full  hits  that  strike  a  particular  user,  we  have 
partial  ones,  as  well;  this  increeises  the  probability  of  a  hit  to  P/,  =  1  —  (1  —  2fq)^~^ .  as 
stated  at  the  beginning  of  this  section.  Substituting  this  for  P/,  in  the  above  expressions 
we  obtain  the  desired  expressions  for  the  multireception  probabilities  P(/,n7  —  l\k). 


3.  A  NEW  APPROXIMATION  FOR  P(hm  -  l\k) 

In  the  previous  section,  we  derived  exact  expressions  for  the  probabilities  P{l.Tn  — 
l\k)  (1  =  0, 1,  •  •  •  ,m,  ni  <  k).  Here  we  develop  an  approximation  method  based  on  the 
Gaussian  multivariate  distribution.  First  we  present  a  general  approximation  technique  for 
P(/,m  —  l\k)  for  systems  with  general  interference  covariance  matrices.  Then  we  exploit 
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the  specific  form  of  the  covariance  matrix  of  the  interference  for  the  symmetric  FH/SSMA 
problem  to  obtain  a  closed  form  expression.  Finally,  we  derive  the  necessary  covariance 
matrices  for  the  various  cases  of  interest. 

3.1  Gaussian  Approximation 

The  computation  of  P{l,m  —  /|^')  is  essentially  a  combinatoric  problem,  which  in¬ 
volves  keeping  track  of  the  bit  erasures  or  errors  and  declares  an  error  when  the  count 
for  a  particular  user  exceeds  the  correcting  capability  of  the  code.  This  is  essentially  the 
multi-dimensional  extension  of  the  binomial  counting  experiment,  which,  however,  results 
in  excessive  computational  complexity,  as  shown  in  the  previous  section.  Here  we  we  ap¬ 
proximate  the  required  “multi-nomial”  probability  distribution  with  a  multidimensional 
Gaussian. 

Let  us  define  the  random  variables  X(,  {\  <  (  <  n,),  where  n  is  the  number  of  bits  per 
packet  for  the  zth  among  m  users  of  interest  (m  <  k)  so  that 

Hi  =  1,  if  the  fth  bit  (or  symbol)  of  user  i  is  incorrect  (due  to  an  erasure  or  an 
error,  depending  on  the  decoding  which  occurs  with  known  probability  p) 
xti  =  0,  if  the  £th  bit  of  user  i  is  correct  with  probability  (1  —  p). 

The  actual  calculation  of  p  presents  no  difficulty  and  will  be  carried  out  in  Section  3.3  for 
the  various  case  of  interest.  Now  define  for  user  i  (1  <  i  <  ni)  the  RV  x,  such  that 

n 

X,  =  ^  xp. 

1=1 

I,  is  the  number  of  bit  errors  or  ereisures  that  sender  i  suffers  among  his  n  bits  within  the 
packet  (slot).  Thus,  0  <  x,  <  n.  xa  is  independent  of  Xft  when  (  ^  F,  because  of  the 
reindom  FH  patterns  assumed.  Consequently,  x<  is  the  sum  of  n  i.i.d.  random  x'ariables 
and  this  tends  towards  a  Gaussian  distribution  for  large  n  [mean  np,  \-ariance  np(l  —  p)]. 
If  we  consider  any  linear  combination  of  the  x,s,  say 


1=1 


a.x, 
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then 


1=1  \/=l 

n  k 

=  X/  X/ 

(=\  i=l 


-E 


It  turns  out  that  the  c^s  are  either  independent  of  each  other  or  one-dependent  (i.e.,  Z(-\, 
Z(,  Z(+i  are  dependent  but  Z(  is  independent  of  2f-|.25 -t-2, 2f-3,  ■  •  Independence 
arises  if  the  system  is  hop  synchronous  and  one-dependence  when  relative  hop  offsets  are 
permitted,  as  in  the  hop-asynchronous  systems  (see  Section  3.3.2).  In  either  case,  7  is  a 
sum  of  i.i.d.  or  one-dependent  RVs  and,  as  n  — »  oo,  tends  to  have  a  Gaussian  distribution. 
Consequently,  all  i,s  are  jointly  Gaussian  if  the  codeword  length  n  is  sufficiently  large. 
Define  the  m-dimensional  column  vectors 


and  //  = 


Then  we  have  the  multiv'ariate  Gaussian  probability  density  function  (pdf) 

p,(x)  =  ,  *-  (3. 

where  S  is  the  m  x  m  covariance  matrix  with  diagonal  elements 

a  =  E{(i,  — (3. 

and  off  diagonal  elements  (all  of  which  are  equal  due  to  the  symmetry  of  our  problem) 

b=  £■  {(i,  -  np)(xj  -  np)}  (3. 


{  ”  IF” 

.^=1 
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a  and  b  are  calculated  in  Section  3.3  on  the  b2isis  of  the  signaling  scheme  (FH  hop- 
synchronous  or  hop-asynchronous)  or  on  the  basis  of  the  decoding  scheme  and  the  presence 
of  AWGN.  We  know  that,  for  user  i  to  be  successful, 


0  <  I,  <  e 

where  e  is  the  erasure  or  error  correcting  capability  of  the  RS  code  (e  =  n  —  kc  or  e  = 
(n  —  kc)l2,  respectively,  for  RS). 

If  user  i  is  unsuccessful, 

e  <  Xi  <  n. 


Hence, 


PG{l,m-\k)=  f  ■■■  f  [  ■  I  p^[x)dxi  ...dxk. 

Jq  Jq  Je+1  Je+1 


In  the  above  equation,  there  are  1  integrals  of  the  form  and  m  —  1  integrals  of  the  form 
f”  We  define 


r  0  1 


and  have 


Pcilm  -  l/k) 


-i: 


^(27r)"'|i:| 


This  integred  can  be  simplified,  if  the  exponent  is  converted  from  a  quadratic  to  a  sum-of- 
squares  form.  This  can  be  materialized  through  a  linear  transformation  that  diagonalizes 
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Now  E  teikes  the  form 


E  = 


abb 

b  a 


b  b  .  a] 


Therefore, 


V  _ 


a  —  b 


a  —  b  0 


a  —  b 


+ 


b  b  . 
b 


bl 


b  ....  bj 


mxm 


mxm 


and,  cons<.*quently. 


^  =  {a- b)I  +  b 


[1  1  1 


mxl 


or,  equivalently, 
where  =  \A[1  1 


E  —  A  -i- 

■ .  l]ixm  and  A=  (a  -  b)I.  Moreover, 

V-.  -  1-1  _ 

'  1  +  v^A~^u 


1 


■I  ~ 


u  u 


{a-b)  (a  -  b)(a (m  -  \)b) 


so  that,  if  we  define 


a  = 


b  = 


a  —  b 


(a  -  b)[a  +  (m  -  l)b] 


(3.6) 

(3.7) 
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then 


b  b  ...  b 
b  b 
b 

b  .  .  b 

One  way  to  diagonalize  E"’  is  to  find  its  eigenvalues  and  the  corresponding  eigenvec¬ 
tors  and  then  create  a  transformation  matrix  with  the  eigenvectors  as  its  columns.  The 
eigenvalues  of  are  given  by  the  equation 

det(A/-E-')  =  0. 


16 


where 


= 


a  —  b 


and 


Ai  = 


1 


a  +  (m  —  l)b 


Finally,  notice  that 


det{T,)  =  (a  —  6)*"  ^  [a  +  (m  —  1)6] 


(3.11) 


3.2  The  Integral  Transform  Method 

Consider  the  exponent  of  the  integrand  in  (3.5).  After  a  shift  of  variables  to  account 
for  the  mean,  it  can  be  simplified  to 


1 


T  —  — 


2{a-b) 


j=i 


a  +  (m  -  1)6  \ 

I  J”"  * 


Using  Uj  =  makes  the  exponent 


1 

0 


J=1 


The  cross  product  terms  come  from  the  factor 

(5- 

But  the  square  in  the  above  equation  can  be  eliminated  with  the  help  of  the  following 
integral  tr2insform  (see  [5]): 


ej"  = 


v/2 


/■+00  yi 

I  dy  e'^^e  2^ 
J  —00 


(3.12) 


TTcr^  J-OQ 


with 


a'  = 


a  -f  (m  —  1)6 


(3.13) 


17 


and 


(f)  =  (3.14) 

J=1 

After  some  manipulations  in  which  we  use  (3.11)-(3.14),  we  obtain  the  basic  result 


?,m  -  l\k)  =  j 


(3.15) 


where 


e  —  np  —  yx/fc 


~np  -  yyjb 
\/a  —  b 


(3.16a) 

(3.166) 


and  e  —  np  is  the  erasure  (or  error)  correcting  capability  of  the  block  code  employed,  shifted 
by  the  meaai  np.  Eq.  (3.15)  gives  a  method  for  calculating  Pa{l,m  —  l\k)  with  linear 
computational  complexity  in  m.  Note  that,  for  the  CEise  6  <  0  (6  being  the  off-diagonal 
elements  of  the  covariance  matrix),  (3.15)  involves  $  evaluated  at  complex  arguments. 
This  is  perfectly  legitimate  and  results  in  an  additional  smcdl  computational  effort.  In  this 
case,  the  real  part  of  the  entity  [$(ri )  —  4>(r2)]^  •  [1  —  ^(-i )  +  ^(■^2)]"*”^  is  involved  in  the 
integral  of  (3.15). 


3.3.  Derivation  of  Received  Code  Vector  Means  and  Covariances 

The  method  developed  in  Sections  3.1  and  3.2  for  approximating  the  probability 
P(/,  m  —  l\k)  requires  only  knowledge  of  the  three  quantities  np,  a,  and  6.  These  quantities 
are  the  meein  of  the  received  code  vector  x  and  the  diagonal  and  off-diagonal  terms  of 
its  covjiriance  matrix  E  [see  (3.0)-(3.3)].  Here  we  evaluate  these  quantities  for  FH  sig¬ 
naling  and  the  various  types  of  system  conditions  (hop-synchronous  or  hop-asynchronous, 
presence  of  AWGN)  and  decoding  schemes  (ereisures-only,  errors-only,  and  errors /erasures 
decoding)  of  RS  codes. 
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Each  transmitter  sends  data  in  packets  with  n  RS  symbols  per  packet.  For  the  jth 
symbol  of  the  zth  user,  the  RV  xji  takes  V7».lue  1,  if  that  symbol  is  either  erased  or  in  error, 
ajid  0  otherwise.  The  probability  of  Xji  being  1  is  p  (to  be  evaluated  separately  later  for 
the  two  cases).  The  expression  Xi  =  counts  the  number  of  errors  or  erasures  per 

packet  and  has  mean  np;  we  also  need  a,  ^ls  defined  in  (3.2)  and  b  as  defined  in  (3.3). 

Each  transmitter  hops  randomly  between  q  available  frequencies  with  one  M-ary  sym¬ 
bol  transmitted  per  frequency  hop.  A  “hit”  takes  place  when  the  frequency  at  which  a 
particular  user  chooses  to  transmit  is  also  used  by  one  or  more  other  users  during  a  dwell 
time.  The  occurrence  of  this  event  can  be  detected  by  listening  to  the  channel  and  con¬ 
ducting  a  threshold  test;  then  a  symbol  erasure  is  declared.  Between  different  users,  the 
duration  of  a  hop  (or  dwell  time)  may  be  assumed  perfectly  synchronized,  or  more  real¬ 
istically,  as  involving  relative  delays.  Noise  (AWGN)  may  also  be  present  in  general,  but 
this  is  dealt  with  a  little  later. 

We  first  evaluate  the  desired  quantities  for  the  noiseless  hop- synchronous  case  in  Sec¬ 
tion  3.3.1,  then  for  the  noiseless  hop- synchronous  case  in  Section  3.3.2;  for  these  cases 
errors-only  zmd  erasures-only  decoding  are  treated  together;  finally,  we  incorporate  the 
effects  of  AWGN  and  treat  errors/ erasures  decoding  in  Section  3.3.3. 

3.3.1  Noiseless  Hop-Synchronous  FH/SSMA 

Since  we  only  have  multiple  access  interference,  the  decoding  employed  is  erasures 
only;  however  the  results  are  also  applicable  to  errors-only  decoding  after  a  trivial  change 
in  the  parameter  e.  The  probability  of  a  hit  is 


and 

a  =  ^{(x, -np)^}  =  np(l-p)  (3.18) 

b=  E{{x,  -np){xj  -np)} 

=  '^{xu-p)  ^(aro-P)  I 
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But 


=  xtj)  -  p^)  • 

e=i 


E{xiixej}  =  Pr{x(,  =  1,  X(j  =  1} 

=  Pr{A^,j,  Xf,  =  1,  X(j  =  1}  +  Pr{Al^,  xe,  =  1,  xij  =  1} 

where  is  the  event  that  users  i  and  j  hop  to  the  same  frequency  during  the  £th  hop. 
Therefore, 


E  {i/.x/j}  =  1  •  Pr{A^}  + 


Y 


Lm=  1 


m  users  hit  i. 


user  j  gets  hit 


1 

9-1. 


k-2 


-f 


(3.19) 

and 

b  =  n  E{x(,  X(j}  —  np^.  (3.20) 

3.3.2  Noiseless  Hop-Asynchronous  FH/SSMA  Channel 

In  Section  3.3.1,  it  was  assumed  that  hops  between  users  were  perfectly  synchronized, 
i.e.,  that  there  was  no  overlap  between  the  Ith  hop  of  user  i  and  the  (/  +  l)th  hop  of 
user  j,  for  any  i,  j.  This  assumption  is  not  realistic,  since  we  assume  that  the  system  is 
distributed  and  not  centredly  controlled.  It  is  time  now  to  relax  that  assumption. 

Assume  that  each  hop  is  of  total  duration  T.  Let  us  set  a  reference  point  for  the 
beginning  of  each  hop.  Then  we  assume  that  the  the  displacement  in  time  of  the  beginning 
of  the  hop  of  any  user  is  eui  RV  uniformly  distributed  between  —Tj2  and  -\-Tj2. 
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First  we  note  that  p,  the  probability  that  any  user  is  hit  while  in  the  Ith  hop,  changes 
as  follows: 


P  = 


(3.21) 


Note  that  m  is  independent  of  Xmt,  for  I  ^  m,  due  to  the  random  assumption  for  the 
hopping  pattern.  However,  because  of  chip  overlap,  xij  is  not  independent  of  imi-  when 
m  cind  /  differ  only  by  1.  Therefore,  we  have 


a  =  E{ix,-np)‘^}  =  np(l-p) 


(3.22) 


b=  E{{x,  -  np){xj  -  np)} 


=  E 


-P) 

.f=i 


-  P) 


=  E{{Xu  ~  Pll-Tij  ~P)  +  (Xu  -  p)(^2j  -  p) 

+  (^2i  -p}  +  (X2,  -pKa-?;  -p)  +  (2'2.  -  P)(2-3j  -  p) 


+  (^^n,  -  p)(2-(„_l);  -  p)  +  (Xni  “  p)[lr,}  -  p)}- 

which,  since  —  p)(x(/+m)j  ~  p)}  =  0  for  any  i,  j  I,  m  =  2,3,....  reduces  to 

b  =  n(pi-p^)  +  (2n -2)(p2 -p')  (3.23) 

where  pi  =  and  p2  = 

Next  we  derive  expressions  for  p\  and  p2.  In  the  following,  whenever  we  refer  to  a  hit. 
we  mean  on  the  /th  hops  of  users  i  and  j. 

Pi  =  Pr{xi,  =  l,X(j  =  1) 

=  Pr{x(,  =  l.Xfj  =  l|i,  j  hit  each  other  on  fth  hop) 

•  Pr{i,  j  hit  each  other  on  fth  hop} 

+  Pr{  xit  =  =  l|i,  j  do  not  hit  each  other  on  Tth  hop} 

■  Pr{i,  j  do  not  hit  each  other  on  f'th  hop} 
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The  first  term  equals  l/q.  The  second  term  is  equal  to 


Pr{x(,  =  =  l|i,  j  do  not  hit  each  other}. 

We  write  the  la>,ter  probability  as 

i  hit  by  G,  i  not  hit  by  6’,  j  is  hit|  i,j  do  not  hit  each  other) 

G 

=  Pr{  i  hit  by  G|  i,  j  do  not  hit  each  other) 

G 

Pr{  j  gets  hit  |  i  hit  by  G.  i.j  do  not  hit  each  other) 

where  G  refers  to  all  possible  groups  of  the  remaining  k  —  2  users  and  all  nits  pertain  to 
the  f'lh  hop. 

Consider  a  particular  user  u  from  group  G  and  the  probability  of  u  hitting  j .  given 
that  u  hits  i  and  i  and  j  do  not  hit  each  other.  We  note  that  there  will  be  two  consecutive 
hops  (dwell  times)  of  it  that  overlap  in  time  with  the  /th  hop  of  i.  Out  of  these  two,  one 
or  both  hit  i.  Therefore,  there  are  two  possible  situations;  (1)  both  hops  of  u  overlap  with 
the  /th  hop  of  j  with  probability  1/2  (given  the  uniform  distribution  of  the  chip  delay  of 
the  users)  and  (2)  only  one  of  these  two  hops  overlaps  with  j  with  probability  1/2.  Note 
that  one  of  the  above  possibilities  has  to  take  place,  since  the  /th  hops  of  ?  and  j  do  overlap 
to  some  extent,  under  the  asynchronous  FH/SS  system  assumption. 

Consider  the  probability  that  a  hits  the  /th  hop  of  given  that  (1)  above  holds,  a 
hits  i,  and  i  and  j  do  not  hit  each  other.  This  probability  is 

11,11  1 

2  '  q  -  1  ^  2  '  q-\  ~  q~\ 


Also,  consider  the  proljabilit}’  that  a  hits  the  /th  hop  of  given  that  (21  above  holds,  that 
a  hits  z.  cUid  that  t  and  j  do  not  hit  each  other.  This  will  be 


111  _  3  1 

~  2’  q-\ 
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From  all  the  above  we  surmise  that  the  probability  of  u  hitting  the  Ith  hop  of  j,  given 
that  u  hits  the  lih.  hop  of  i  and  that  i  and  j  do  not  hit  each  other  on  the  /th  hop,  is 

11  13  1  5  1 

2'  q~  I  2’  2'  q -I  ~  4  q-l' 

Finally, 

Oi  =  Pr{xi,  =  l,xij  =  ll  i  j  do  not  hit  each  other} 


(3.24) 

and  thus 

px  = 

We  now  derive  p2  ~  Here,  when  talking  of  a  hit.  we  refer  to  the  1th 

hop  of  i  and  to  the  (1  +  l)th  hop  of  j.  We  have 

P2  =  Pr{xu  =  1  X(,+  i)_,  =  1} 

=  Pr{x/,  =  1,X(;4.2)^  =  1|  i,  j  overlap  in  time) Pr{  i,j  overlap  in  time} 

-f  Pr{xi,  =  l,X(/^.i)j  =  Ij  i,  j  do  not  overlap  in  time}Pr{  i,  j  do  not  overlap  in  time} 

In  the  above  equation, 

Pr{  i,j  overlap  in  time}  =  Pr{  i,  j  do  not  overlap  in  time}  =  1/2 

and 

Pr{x(,  =  l,X(;+i)j  =  Ij  i,  j  do  not  overlap  in  time}  =  Pr{x(,  =  l}Pr{x(;+,)j  =  1}  = 


23 


Moreovei , 


Pr{xi,  =  =  1|  i,  j  overlap  in  time} 

=  P7'{xi,  =  l,X(^i^i)j  =  1|  i)  j.  overlap,  hit  each  other  on  f  +  1th  hop) 

•  Pr{  i,j  hit  each  other  on  ^  +  1th  hop|i,  j  overlap} 

+  Pr{xi,  =  l,X(j^.i)j  =  Ij  i,  j  overlap,  do  not  hit  each  other  on  ^  +  1th  hop} 

•  Pr{  i,  j  do  not  hit  each  other  on  ^  +  1th  hopji,  j  overlap} 

From  this  point  on,  the  evaluation  of  the  quantities  of  interest  proceeds  in  a  similar  way 
to  that  for  pi  and  thus  it  will  be  omitted.  It  turns  out  that 

02  =  P7'{xi,  ~  =  1|  i,  j  overlap,  do  not  hit  each  other}  =  oj  (3.26) 

where  oi  is  given  by  (3.24),  so  that  the  final  result  is 

P'2  =  \  -  |  +  (3.27) 

The  desired  quantity  b  is  now  obtained  from  (3.24)-(3.25)  and  (3.26)-(3.27)  upon 
substitution  in  (3.23). 

3.3.3  Inclusion  of  AWGN  at  Front  End  of  Each  Receiver 

So  far,  we  have  assumed  that  our  channel  is  noiseless  and  that  we  can  employ  erasures- 
only  decoding,  sir^c  symbol  erasures  can  be  detected.  If  we  incorporate  noise,  the  possi¬ 
bility  of  a  symbol  error  exists.  Symbol  error  can  only  take  place,  if  there  is  no  erasure  due 
to  multiple-access  interference  or  to  hits.  Now  we  shall  employ  errors/erasures  decoding. 
We  maintain  for  user  the  RV  c/,  which  is  set  to  1,  if  the  Ith  bit  is  an  erasure,  and  the 
wi,,  which  is  set  to  1,  if  it  is  in  error.  Define  i/,  =  c/,  -|-  2u’/,.  Let 


We  know  that,  for  user  ?  to  be  succes.sful. 

0  <  X,  <  u  —  kc 
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Using  arguments  similar  to  those  of  the  noiseless  case,  we  can  prove  once  again  that  the  x,s 
are  jointly  Gaussian,  each  being  the  sum  of  one-dependent  variables.  We  need,  therefore, 
to  find  £'{ij};  foi  this  we  need  E{wit}  zmd  the  probability  of  M-ary  symbol  error  ps, 
which  is 


PN  =  Pn{M) 


E 

m=l 


M-l 

m 


+  l  ,n  rlog^MEt 
■  g  nv^T  A'  Q 

m  -f  1 


(3.28) 


where  r  is  the  code  rate  and  Ei,/Nq  the  signal-to-noise  ratio.  Then 

£{r.}  = 

=  n  •  [p  -I-  2  •  p.v  •  (1  -  p)j 

=  np  (3.29a) 


where 

p'=p  +  (l-p)p^  (3.296) 

and  p  is  given  by  (3.21)  or  (3.17)  for  hop-asynchronous  or  hop-synchronous  FH/SSMA 
systems,  respectively. 

Note  that  the  probability  of  an  erasure  E{zi,}  does  not  change  with  the  presence  of 
noise. 

Further  we  need  to  calculate  a'  which  is 
a'  =  £{x,2}-£{x,}" 

=  n  ■  ^E{{zi, +2wi,)^}  ~  p'^ 

,2 


(3.30) 


—  n-  E{zi,^  +  4  u’h^}  +  4  E{zi,wi,} — 


=  n  ■ 


P+  4  (1  -  p)  •  p/v  -  p 


We  also  need  6^  which  is  given  by 


6'  =  n(p/-p'  )  -I-  (2n  -  2)(p2' -p'^)  (3.31) 

and  for  whose  computation  p/  and  p'/  ar**  necessary.  The  fo^mt  i  is  defined  as 
Pi'  =  E{xt,xij)  =  E{zi,zij  4- 4zi,wij  +  4wi,wij}. 
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E{ziiZij}  is  the  same  with  that  of  the  noiseless  case,  which  was  dealt  with  in  the  previous 
section  and  found  to  be  equal  to  pi ,  which  is  cedculated  there.  Moreover, 


E{ziiWij}  =  Pr{zu  =  =  1} 

=  Pr{zu  =  1,2/j  =  0}  •  9 
=  [Pr{zi,  =  1}  -  Pr{zii  =  l,zij  =  1}]  Pn 
=  {P-  Pl)PN 

where  pi  is  given  by  (3.25).  Finally, 

E{wi,wij}  =  Pr{wu  =  l.uuj  =  1} 

=  Pr{zn  =  0,2ij  =  0} 


Combining  all  the  above  we  obtain 


/ 

Pi 


pi  +  4(p  -  Pi  )pA'  +  4 


(3.32) 


The  other  quantity  of  interest,  p2*  =  is  calculated  in  an  identical  manner 

and  is  found  to  be 


P2'  =  +  4u>;,u;(/+i)j) 


where 

E{zi,Z^i^l)j  =  P2 

=  Pr{z{,  =  =  1} 

=  Pr{zt,  =  l,2(j+i)j  =  0}  ■  9 
=  [Pr{zi,  =  1}  -  Pr{zi,  =  1,  Z(i+i)j  =  1}  •  PA 
=  (P  -  P2)PN 


E{wiiW[i+i^j}  =  Pr{wii  =  =  1} 

=  Pr{zii  =  0,2(,+i)j  =  0} 
=  Pr{wn  =  l,wij  =  1} 


Combining  the  above  we  obtain 


P2 


p2  +  4(p  —  ^2)PN  +  4 


fc-2 


p\- 


(3.33) 


The  desired  quantity  h'  is  now  obtained  from  (3.32)-(3.33)  and  (3.29)  upon  substitution  in 
(3.31).  and 

The  above  results  are  valid  for  erasures/errors  decoding,  when  information  about 
the  state  of  the  channel  (side  information)  pertaining  to  hits  from  other  users  is  available. 
If  error-only  decoding  is  employed  (in  the  absence  of  side  information),  the  results  are 
simplified  considerably.  In  that  case,  we  use  the  previous  definition  of  x/,  to  denote  the 
presence  of  a  error  (not  an  erasure)  at  the  /th  hop.  Then  the  necessary  quantities  p',  o', 
and  b'  are  obtained  from  the  quantities  p,  pi,  and  p2  of  the  asynchronous  case  (no  AWGN) 
of  the  previous  section  as 

p=p  +  (l-p)pN,  (3.34) 

where  p  is  given  by  (3.21)  or  (3.17)  and  pN  is  given  by  (3.28), 


a  =  np'(l  -  p'). 


(3.35) 


and 

b'  =  n[pi  +  (1  -  P\)pn  ~P^\  +  (2n  -  2)[p2  +  (1  -  P2)Pn  -  P^]-  (3-36) 


4.  APPROXIMATION  BASED  ON  THE  INDEPENDENCE  ASSUMPTION 

The  eissumption  of  independence  between  the  packet  errors  of  the  users  is  commonly 
made  for  simplifying  the  evaluation  of  P{l,m  —  l\k).  The  relevant  expressions  are  given  here 


for  reference,  because  they  will  be  used  for  computing  the  numerical  results.  Let  e  and  t  be 
the  erasure  and  error  correcting  capability,  respectively,  of  the  block  code  employed,  e  the 
probability  of  symbol  erasure,  p'  the  probability  of  symbol  error  for  errors-only  decoding, 
p  p^cbabilit;’  cf  symbol  error  for  errors  and  erasures  decoding,  q{M)  the  probability  of 
an  M-ary  FSK  symbol  error  due  to  AWGN,  and  Pe(^)  the  probability  of  a  codeword  (or 
receiver)  error  for  the  typical  user.  Also,  let  P„d{l,  rn  —  be  the  probability  that  exactly 
m  —  I  users  suffer  packet  error,  given  that  k  users  transmit,  and  that  the  probability  is 
computed  using  the  independence  assumption.  Then  we  have 

P^^d^T[Um-l\k)  ^  {^^P.rAUm-l\k)  (4.1) 

and 

P„d{Um-l\k)=  \p,{k)r-‘  [l-p,{k)]‘.  (4.2) 

where 

P.(t)  =  y  (")('■  (!-<)"-'■  (4.3) 

l.iTu 

with  e  —  n  —  kc  for  erasures-only  decoding; 

=  E  (,")/'  (!-?')"■''  (4.4) 

l=.+i 

where  t  =  (n  —  fcc)/2,  for  errors-only  decoding;  emd 

n  n  —  li 

p.(i-)  =  E  E 

/i=0  f2=:max{e+l-2/i,0) 

where  e  =  n  —  k^  for  errors  and  erasures  decoding  For  example  for  asynchronous 
FH/SSMA  with  AWGN  the  parameters  c,  p',  and  p  take  the  \’alues 


n  -li 

h 


p-e) 


n-/i  -/j 


e=  1  - 


(4.6) 


p  =t  +  {l  -  e)pN 


(4.7) 


and 


p  =  (1  -  c)Pjv 


(4.8) 


Vkhere  pyv  is  giver  by  (?.13). 

We  can  show  that,  under  light  traffic  conditions  (i.e.,  k/q  <<  1),  the  approximation 
based  on  the  independence  assumption  (4.2)  is  valid  for  all  /,  m,  and  k.  To  this  end  it 
suffices  to  show  that  P(P„„)  =  Ps{p,m  —  p.,k,q)  of  (2.8)  can  be  approximated  by  the 
corresponding  expression,  when  all  receivers  operate  independently.  Indeed,  for  k  <<  q 
we  have 


Ps{p,rn  -  p,9.  A')  ~  [1  -  Ph{q,k)Y 


>=]  ^  '  j-i 


=  [1  -  P^(q,A-)]"  ■  [1  -  (1  -  Ph{q.k)T~‘‘ 

=  [l-P,(9,A-)l^Pfc(g,A:)"‘-"  (4.9) 


5.  NUMERICAL  RESULTS 

In  this  section  we  present  our  numerical  results  and  comparir-ous.  As  explained  before, 
the  exact  numerical  evaluation  of  P(/,  m  —  l\k)  becomes  prohibitive  for  m  >  4  and  so  only 
results  for  m  =  2  and  m  =  3  can  be  presented.  Indeed,  as  it  is  clear  from  equation  (3), 
2^"  —  1  nested  sums  are  necessary  for  the  computation  of  P{l,m  —  l\k)  in  the  case  of  errors- 
only  or  erasures-only  decoding,  whereas  3”*  —  1  such  sums  are  necessary  for  the  case  of 
errors/ erasures  decoding  [see  (17)].  Since  in  each  sum  the  number  of  terms  v’aries  between 
0  and  #  or  #  +  1  and  n,  zis  many  2is  (t  -j- 1)^"* or  even  (#  +  1)^  terms  may  be  necessary. 
Eissuming  n  =  2t.  This  implies  that,  for  a  (32, 16)  RS  extended  code  and  m  =  2  receivers. 
3  and  8  nested  sums  are  necessary,  respectively,  for  a  toted  of  17^  or  17®  terms;  while  for 
m  =  3,  7  and  26  nested  sums  eire  necessary  for  a  maximum  of  17'  or  17^®  terms.  The 
corresponding  numbers  for  m  =  4  are  15  and  80  nested  sums  (i.e.,  a  maximum  of  17'^  and 
17®°  terms),  which  eire  prohibitive  for  most  computers. 

All  results  presented  in  this  section  pertain  to  eisynchronous  FH/SSMA  systems  em¬ 
ploying  RS  (32, 16)  extended  codes  with  32-ary  FSK  data  modulation  and  noncoherent 
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demodulation.  Thus  M  =  n  =  32  and  kc  =  16  for  the  modulation  and  code  parameters  of 
interest.  The  performance  of  different  minimum-distance  decoding  algorithms  is  evaluated 
including  errors-only,  erasures-only,  cuid  errors/erasures  correction  decoding.  The  effect  of 
AWGN  is  taken  into  account  in  all  cases.  The  number  of  frequencies  available  for  hopping 
is  q  and  m  is  the  number  of  FH/SS  receivers  of  interest.  In  all  cases,  results  based  on  the  ex¬ 
act  expressions  P{l,'m  —  l\k)  and  results  based  on  the  independence  approximation  (IROA) 
P/(/,m  —  l\k),  or  the  Gaussian  approximation  PG(/,m  —  l\k),  are  presented.  Moreover,  the 
results  for  the  packet  (codeword)  probabilities  P(/,  m  —  /|fc)  and  P/(/,  m  —  l\k)  are  presented 
in  different  subtables  as  the  number  of  correct  packets  I  changes  (Z  =  m,  m  —  1,  •  •  • ,  1, 0). 

In  Tables  1  and  2  we  compare  the  performance  of  m  =  2  and  m  =  3  receivers  employing 
different  forms  of  minimum-distance  decoding,  in  particular  errors-only,  erasures-only,  and 
errors/erasures  decoding,  for  different  values  of  the  information  bit  signal-to-noise  ratio 
Eb/No-  The  total  number  of  contending  users  is  A:  =  10  and  the  number  of  frequencies 
q  =  100.  The  superiority  of  erasures  decoding  and  errors/erasures  decoding  over  errors- 
only  decoding  is  established  for  the  range  of  values  of  Eb/No  considered.  The  limiting 
values  of  the  packet  probabilities  are  essentially  achieved  already  for  Eb/No  =  10  dB,  at 
which  point  all  errors  are  caused  by  other-user  interference.  The  approximation  based  on 
the  IROA  appears  to  be  very  close  to  the  exact  result  for  most  cases.  The  accuracy  of  the 
IROA  appeeirs  to  be  better  for  P(2, 0|A')  and  P(3, 0|A.')  than  it  is  for  P(l,  1|A')»  -^(2, 1|A:)  or 
P(l,  2|A:);  it  is  also  better  for  errors-only  or  erasures-only  decoding  than  for  errors/erasures 
decoding.  Finally,  accuracy  improves,  as  Eb/No  increases.  The  exact  values  of  P{l,m  —  l\k) 
are  missing  from  Table  2  because  the  computation  is  prohibitive;  for  m  =  3,  26  nested 
sums  and  a  maximum  of  17^®  terms  are  necessary  in  (17). 

Tables  3  and  4  illustrate  the  performance  of  rr?  =  2  and  m  =  3  receivers,  respectively, 
for  FH/SSMA  systems  employing  errors-only  decoding,  as  the  total  number  of  users  k 
and  hopping  frequencies  q  vary.  The  information  bit  signal-to-noise  ratio  Eb/No  =  10 
dB.  Comparing  the  exact  results  with  those  obtained  under  the  IROA  we  observe  that 
the  accuracy  of  the  latter  is  better  for  P(2,01A:)  and  P(3,0|A')  than  for  the  rest  of  the 
P{l,m  —  l\k)s.  For  fixed  q,  as  the  number  of  contending  users  k  increases,  the  approximation 
becomes  less  accurate.  However,  for  large  values  of  k,  both  P{l,m  —  /(A:)  and  Pi{l,m  —  /|A-) 
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become  very  small.  On  the  other  hand,  for  fixed  k,  the  accuracy  improves  as  q  increases. 
Overall,  the  accuracy  of  the  IROA  improves  as  the  ratio  k/q  decreases.  This  verifies  the 
analysis  of  Section  3  under  light  traffic  conditions. 

Similarly,  Tables  5  euid  6  illustrate  the  performeince  of  m  =  2  and  m  =  3  receivers, 
respectively,  for  FH/SSMA  systems  employing  erasures-only  decoding,  as  the  total  num¬ 
ber  of  users  k  and  hopping  frequencies  q  veiry.  The  information  bit  signal-to-noise  ratio 
Eb/No  =  10  dB.  Similar  observations  to  the  ones  made  for  Tables  3  and  4  are  valid  here. 
Moreover,  the  results  of  Tables  5  and  6  are  uniformly  better  than  those  of  Tables  3  and 
4,  since  erasures-correction  decoding  is  considerably  more  powerful  than  errors- correction 
decoding  for  for  RS  codes  (RS  codes  can  correct  twice  as  many  erasures  as  errors). 

An  interesting  fact  that  holds  true  in  the  numerical  analysis  we  have  performed  is 
that  P{m,0\k)  is  larger  than  P7{m,  0|t)  in  all  cases.  In  other  words,  IROA  seems  to  give 
“optimistic”  results  in  comparison  to  the  exact  analysis. 

In  Tables  7  to  10  the  exact  results  for  P{l,m  —  /|/r)  are  compared  to  those  obtained 
via  the  Gaussian  approximation  method  described  in  Section  3,  the  latter  is  denoted  by 
Pc{Um  —  l\k).  The  same  system  assumptions  presented  at  the  beginning  of  this  section 
hold.  As  before  the  results  for  the  packet  (codeword)  probabilities  P(/,m  —  l\k)  and 
PG{l,m  —  l\k)  are  presented  in  different  subtables  as  the  number  of  correct  packets  I 
changes.  There  is  actually  a  one-to-one  correspondence  between  Tables  1  and  7  and  Tables 
2  and  8.  There  are  also  similarities  in  the  organization  of  the  results  between  Tables  3 
and  9  and  Tables  4  and  10.  However,  in  Tables  9  and  10  the  number  of  simultaneous 
transmission'-  k  changes  whereas  the  number  of  frequencies  is  held  constant  at  5  =  100 
and  the  results  axe  presented  for  all  three  types  of  decoding  considered  in  this  report: 
errors-only,  errors/erasures,  and  erasures-only  decoding. 

As  it  becomes  clear  from  Tables  7  to  10  the  Gaussian  approximation  is  not  as  close 
to  the  exact  results  as  the  IROA  is,  at  least  for  the  cases  of  m  =  2  and  m  =  3  receivers 
and  rj  =  32.  It  appears  that  the  accuracy  of  the  approximation  improves  as  the  number 
of  simultaneous  users  k  increases  aiad  in  general  it  is  better  for  larger  values  of  the  exact 
probabilities. 
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6.  CONCLUSIONS 

For  FH/SSMA  communications,  we  presented  exact  expressions  for  the  multireception 
probabilities  P{1,  m  —  /|fc)  for  hop-synchronous  and  hop- asynchronous  systems  and  various 
types  of  RS  decoding.  The  effects  of  AWGN  were  taken  into  consideration.  However,  these 
expressions  are  very  difficult  to  evaluate  for  m  >  4,  as  they  require  computation  of  2"’  —  1 
or  3"*  —  1  sums. 

We  also  established  the  validity  of  the  IROA  for  the  case  q  >>  k.  Additionally,  our 
numerical  analysis  indicated  that  IROA  is  a  good  approximation  for  the  multireception 
probabilities,  for  m  —  2  and  m  ~  3.  The  accuracy  of  the  approximation  depends  on  the 
specific  values  of  q  and  k]  it  improves  with  decreasing  k/q;  for  q  >>  k,  the  corresponding 
results  are  almost  identical  to  the  exact  ones.  Therefore,  it  appears  that  IROA  gives  useful 
results,  while  requiring  minimal  numerical  effort. 

The  kind  of  behavior  observed  so  far  is  expected  to  be  similar  for  higher  values  of  m,  as 
well.  However,  this  is  only  a  conjecture  at  this  time  amd  additional  work  on  the  derivation 
of  computationally  efficient  techniques  for  the  evaluation  of  P(/,m  —  l\k)  is  necessary  to 
prove  this  claim. 

Note  that,  in  many  practical  applications,  the  generated  traffic  is  light.  In  such  cases, 
the  condition  q  >>  k  is  eaisily  satisfied  and  the  IROA  can  be  used  to  obtain  the  multirecep¬ 
tion  probabilities.  However,  there  are  also  many  practical  situations,  in  which  q  and  k  are 
of  comparable  magnitude  or  even  k  >  q  and  in  which  the  IROA  can  not  be  validated.  In 
such  situations,  another  method  for  evaluating  the  multireception  probabilities  is  needed, 
since  the  values  of  m  that  are  of  interest  can  be  considerably  larger  than  3  (the  practical 
computational  limit  of  the  exact  approach  presented  in  this  report). 

In  our  report  we  also  introduced  an  approximation  based  on  central  limit  theorems 
for  multiv2Lriate  distributions.  This  Gaussian  approximation  has  low  computationeJ  com¬ 
plexity  and  promises  to  have  good  accuracy  for  large  n  (number  of  symbols  per  codeword 
or  packet).  Unfortunately,  as  our  extensive  comparisons  with  the  exact  results  indicated 
in  Section  5,  for  the  nominal  value  n  =  32  the  Gaussian  approximation  does  not  yield  uni¬ 
formly  satisfactory  results.  We  conjecture  that  the  accuracy  of  the  approximation  improves 
as  n  increases  but  we  can  not  prove  this  since  the  computational  effort  for  e\'aluating  the 
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exact  results  is  prohibitive  for  large  n. 

Finally,  notice  that  the  approach  presented  in  this  report  for  obtaining  the  exact  results 
is  only  applicable  to  FH/SS  systems.  Both  the  IROA  and  Gaussian  approximations  are 
applicable  to  direct-sequence  (DS)  MA  systems  but  an  approach  for  the  exact  evaluation 
of  P(/,m  —  l\k)  is  needed  in  order  to  validate  the  accuracy  of  the  two  approximations  in 
the  DS  case 
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APPENDIX  A 


Al.  Derivation  of  (2.3) 

Here  we  derive  the  expression  for  P(/,  ni  —  /|A:)  as  a  function  of  the  probabilities  Pe...c- 
Let  q,  k  denote  the  number  of  frequency  slots  and  contending  users  in  the  slot,  respectively. 

For  receiver  z,  1  <  i  <  m,  let  be  a  vector  having  0  in  the  positions  of  correct  symbol 
reception  and  1  in  the  positions  of  erroneous  reception,  that  is. 

.^1  ~  (^il  1  ^i2i  ■  ■  ■  >  e|n). 

For  hard  decision  decoding,  the  itk  receiver  decodes  correctly  the  received  packet  iff 

n 

0  <  t 


while  it  decodes  erroneously  iff 


n 

t  -h  1  <  ^  e,j  <  n. 
j=i 

Let  us  now  turn  our  attention  to  the  interreceiver  operation.  For  the  jth  transmitted 
symbols,  for  1  <  j  <  n,  we  define  the  vector  event  Ej  as 

—  i^lji  ^2ji  ■  ■  ■  1  Cm;)- 

As  each  e,j  takes  on  two  possible  values,  there  is  a  total  of  2"‘  possible  Ej,  for  each 
j.  As  slotted  operation  is  assumed  throughout  the  paper,  statistics  Eire  the  same  from 
symbol  to  symbol,  so  that  the  description  of  the  system  is  independent  of  the  particular 
symbol  j.  Thus  we  can  drop  the  dependence  from  j  in  our  notation  and  arbitrarily  assign 
events  E  to  vectors  of  symbol  events  (e; ,  62,  ■  •  • ,  Cm)-  However,  for  simplicity  we  choose 
the  correspondence 

^nn  —  (ci,e2,’**?Cm) 

so  that  (ei ,  C2i  •  •  • ,  Cm)  is  equal  to  the  binary  representation  of  nn.  If  we  define  by  the 
number  of  times  a  particular  event  occurs  eis  the  index  of  symbols  within  a  codeword 


j  varies  (1  <  j  ^  we  obtain,  due  to  the  memoryless  operation  assumption,  the  result 
given  by  (2.3). 

A2.  Solution  to  Recursion  (2.6) 

Here  we  prove  that  Pj(0,  m,q,  k),  given  in  (2.7),  is  the  solution  to  the  recursive  equa¬ 
tion  described  by  (2.6).  For  compactness,  we  denote  the  binomial  coefficients  f  ^  1 


Cn,,.. 


First  we  observe  that,  for  the  binomial  coefficients  Cm.n  the  following  recursion  is  true 


Cm-t-l,!  —  ^m.i— 1  +  c  m  .1  - 


Then,  by  direct  substitution  of  (  2.7)  to  the  right  hand  side  of  (2.6),  we  get 


1  +  ^(-1)‘C„.-,,  ■nil  -  fhiq  -  j  +  l.k  -  j  -i-  1)] 


-[1-P,(q,^)]n+  ^(-1)'C,„-,,.  ni'  -Ph{q~j-k  ■;)][  = 


=  l  -  m[l  -  Phiq,k)]  +  ^(-1)'  (Cm-I.i  +  Cm-i,i-i)  J^[l  -  P/i(<?  -  J  +  T  +  1)1 


+  (-l)'"C,„_l,,„_,  fj[l  -  P,{q  -J  +  I,k-J  +  1)1 


=  1  +  ^(-1)’C„,,.  fill  -PH[q~J+hk-j  +  1)] 


(A2.2) 


From  (A2.2)  we  see  that  (2.6)  has  as  solution  the  expression  given  in  (2.7). 


.^6 


T 


Table  1 

Error  probabilities  computed  using  exact  and  IROA  models  for 
asynchronous  FH/SSMA  with  RS  (32.16)  coding.  ^=10  and  ^=1(X). 

(a)  Probabilities  P  (2,0 1  it)  (exact)  and  Pj{2X)\k)  (under  IROA) 


errors  decoding 

errors/erasures  decoding 

erasures  decoding 

(dB) 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

6 

0.6320 

0.6255 

0.9900 

0.9900 

9.9998 

0.9998 

8 

0.8490 

0.8465 

0.9999 

0.9999 

0.9999 

0.9999 

10 

0.8633 

0.8610 

0.9999 

0.9999 

0.9999 

0.9999 

12 

0.8634 

0.8612 

0.9999 

0.9999 

0.9999 

0.9999 

JO 

0.8635 

0.8612 

0.9999 

0.9999 

0.9999 

3.9999 

(b)  Probabilities  P(  1.1  I  A  )  (exact)  and  P/(l.l  U)  (under  IROA) 


errors  decoding 

errors/erasures  decoding 

erasures  decoding 

£(,7^0  (dB) 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

6 

0.1588 

0.1653 

4.176jrl0“^ 

4.983x10"^ 

5.777x1  Or^ 

5.778x1 0"^ 

8 

T.IOUIO'^ 

7.354i  10'^ 

3.249x10-^ 

1.055x10'^ 

3.341x10"* 

3.341x10"® 

10 

6.462;tl0"^ 

6.687X  10"2 

2.510x10'^ 

2.558x10"* 

2.511x10"* 

2.51x10"* 

12 

6.454X  10'^ 

6.678;cl0"^ 

2.501x10“^ 

2.502X  10"* 

2.501x10"* 

2.501x’0'* 

oo 

6.454;r]0”^ 

6.678J10"* 

2.501J  10"* 

2.502x10"* 

2.501x10"* 

2.501x10"* 

(c)  Probabilities  P{Q2\k)  (exact)  and  P/(0.2l/()  (under  IROA) 


£fc//^o(dB) 

errors  decoding 

errors/erasures  decoding 

erasures  decoding 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

6 

8 

10 

12 

oo 

5.018x10  ' 
8.915x10"^ 
7.438X  10'^ 
7.421x10"^ 
7.421x10"^ 

4.371x10"'^ 
6.389X  10"^ 
5.193x10"’ 
5.179x10"’ 
5.179x10"’ 

9.901X  10"’ 
1.460x10"’ 
9.637X  10'* 
1.458x10"'° 
7.616x10'" 

2.508xi0"^ 
1.113x10"'° 
6.545x10"'" 
6.258x10"” 
6.252X  10"'" 

1.284x10"* 
1.198x10"'° 
7.660x10"" 
7.616x10“" 
7.61 6x  10"" 

3.339X  10"^ 
1.116x10'" 
6.305X  10"'^ 
6.258x10"'" 
6.258x10"'’ 
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Table  2 


Error  probabiliiies  computed  using  exact  and  IROA  models  for 
asynchronous  FH/SSMA  with  RS  (32,16)  coding,  ^=10  and  <7=100. 

(a)  Probabilities  P(3,0I/:)  (exact)  and  P/(3,0l<t)  (under  IROA) 


E,/No  (dB) 

errors  d 

ecoding 

errors/erasures  decoding 

erasures  decoding 

Exact 

IROA 

IROA 

Exact 

IROA 

6 

0.5099 

0.4947 

0.9851 

0.9998 

0.9998 

8 

0.7857 

0.7788 

0.9999 

0.9999 

0.9999 

10 

0.8051 

0.7990 

0.9999 

0.9999 

0.9999 

12 

0.8054 

0.7993 

0.9999 

0.9999 

0.9999 

oo 

0.8054 

0.7992 

0.9999 

0.9999 

0.9999 

(b)  Probabilities  P(2,l  \k)  (exact)  and  P/(2,l  \k)  (under  IROA) 


J7  IKI 

errors  decoding 

errors/erasures  decoding 

erasures  decoding 

Exact 

IROA 

IROA 

Exact 

IROA 

6 

0.1121 

0.1308 

4.958X  10-^ 

5.776x10-^ 

5.778X  10“* 

8 

6.334x10-2 

6.766x10-2 

1.055x10-^ 

3.340x10“* 

3.341x10“* 

10 

5.816x10-2 

6.205X  10-2 

2.558x10"^ 

2,511x10“* 

2.511x10“* 

12 

5.809X  10-2 

6.198x10-2 

2.50  lx  10-^ 

2.501x10“^ 

2.501X  10“* 

OO 

5.809X  10-2 

6.198x10-2 

2,501x  10-* 

2.501x10“* 

2.501X  10“* 

(c)  Probabilities  P  (1,2  II:)  (exact)  and  P/(1,21I:)  (under  IROA) 


Eh/No  (dB) 

errors  decoding 

errors/erasures  decoding 

erasures 

decoding 

Exact 

IROA 

IROA 

Exact 

IROA 

6 

3.678X  10-2 

3.457X  10-2 

2.496X  10-2 

1.284x10-* 

3.339X  lO-*^ 

8 

7.67x  10-^ 

5.878X  10-2 

1.1 13x  10“‘° 

1.198x10-'° 

1.116x10-“ 

10 

6,466x10-2 

4.819x10-2 

6.545X  10- '2 

7.660x10-“ 

6.305x10“'2 

12 

6.452X  10-^ 

4.806x10-2 

6.258X  10-'2 

7.615x10-“ 

6.258x10“ '2 

OO 

6.452X  10-2 

4.806x10-2 

6.258X  10-'2 

7.615x10“" 

6.258x10“' 2 

(d)  Probabilities  F (0,3  Ik)  (exact)  and  P/ (0,3  Ik)  (under  IROA) 


Eb^No  (dB) 

errors  decoding 

errors/erasures  decoding 

erasures  decoding 

Exact 

IROA 

IROA 

Exact 

IROA 

6 

1.340x10“^ 

9.139x10-2 

1.256X  10“^ 

6.117x10-'2 

1.929x10“'2 

8 

1.242x10-2 

5.107x10-^ 

1. 175x10-'* 

7.264X  10“'2 

3.728x10“‘2 

10 

9.718x10^ 

3.742x10"* 

1.674x10“ '2 

1.169x10-‘2 

1.583xI0-'2 

12 

9.687x10“^ 

3.727x10"' 

1.566X  10-'2 

8.808X  10“'2 

1.565x10“'2 

OO 

9.687x10^ 

3.727x10“^ 

1.566x10“ ‘2 

9.580X  10“'2 

1.565x10“'2 
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Table  3 


Error  probabilities  computed  using  exact  and  IROA  models  for 
asynchronous  FH/SSMA  with  RS  (32,16)  coding  and  error-correction  decoding 

(£i//^o  =  10  dB) 

(a)  Probabilities  P  (2.0 1  it)  (exact)  and  P/(2,0l/t)  (under  IROA) 


k  = 

=  5 

k  = 

10 

k  = 

20 

k  = 

50 

9 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

A  Ax  10~® 

9.6j:10'^ 

1.6;c10-2^ 

1.3x10"^^ 

0.0 

0.0 

0.0 

0.0 

50 

0.9200 

0.9172 

0.1065 

0.0967 

9.07X  10“'' 

7.06X  10'^ 

2.4x  lO-^’ 

1.9x10-2’ 

100 

0.9989 

0.9989 

0.8635 

0.8612 

7.36x  10~2 

6.98x  10-2 

2.6x  10-‘° 

2.3x  10-’° 

(b)  Probabilities  P(l,l  i/k)  (exact)  and  P/(l,l  li)  (under  IROA) 


k  = 

=  5 

k  = 

10 

k  = 

20 

k  = 

50 

<7 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

^8x10-5 

9.8x10"^ 

4.1x10“’^ 

3.7x  10-’^ 

0.0 

0.0 

0.0 

0.0 

50 

3.77x  10-2 

4.05X  10-2 

0.2045 

0.2142 

8.39x10-^ 

8.39x  10-^ 

4.9x  10-'^ 

4.3x  10-’^ 

100 

5.45X  10^ 

5.51x10-^ 

6.45x  10'2 

6.67x10-2 

0.1906 

0.1944 

1.5x10"^ 

1.5x10-* 

(c)  Probabilities  P(0.2  l/t)  (exact)  and  P/(0,2l/c)  (under  IROA) 


k  = 

=  5 

k  = 

10 

k  = 

20 

k  = 

50 

9 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

0.999 

0.999 

0.999 

0.999 

1.0 

1.0 

1.0 

1.0 

50 

4.64 -t  10"* 

1.79x10-* 

0.4846 

0.4748 

0.9983 

0.9983 

0.999 

0.999 

100 

6.09X  lO"^ 

3.04x10-2 

7.42x10-* 

5.18x10"* 

0.5451 

0.5413 

0.999 

0.999 
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Table  4 


Error  probabilities  computed  using  exact  and  IROA  models  for 
asynduoncus  FH/SSMA  willi  RS  (32,16)  coding  ana  error-correction  decoding 

(£fcWo=  10  dB) 

(a)  Probabilities  P(3,0I/:)  (exact)  and  P/ (3,0 lit)  (under  IROA) 


k  = 

=  5 

k  = 

10 

k  = 

=  20 

k  = 

50 

<? 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

9.3x10-“ 

9.4x10-“ 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

50 

0.8863 

0.8784 

3.97x  10-2 

3.0x  10-2 

1.2x10-® 

5.9x  10-’° 

0.0 

0.0 

100 

0.9983 

0.9983 

0.8054 

0.7993 

2.1x10“2 

1.8x10-2 

4.9x  10-“ 

3.4x  10-“ 

(b)  Probabilities  Pil.Wk)  (exact)  and  P/ (2,1  lit)  (under  IROA) 


k  = 

=  5 

k  = 

10 

k  = 

=  20 

k  = 

50 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

4.35x  10-* 

9.62x10-® 

1.6x10-2® 

1.3x10-2® 

0.0 

0.0 

0.0 

0.0 

50 

3.37x  10-2 

3.88x10-2 

6.7x10-2 

6.7x10-2 

9.1x10-2 

7.1x10-2 

2.4x10-2® 

1.8x10-2® 

100 

5.39x  10^ 

5.51x10"' 

5.8x  10-2 

6.2x  10"2 

5.2x  10-2 

5.14x10-2 

2.5x10-'° 

2.2x  10-'° 

(c)  Probabilities  P(h2\k)  (exact)  and  P/(l,2 1^)  (under  IROA) 


k  = 

=  5 

k  = 

10 

k  = 

20 

k  = 

50 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

9.8x  10-^ 

9.8x  10-2 

4.  lx  10-'2 

3.7x10"’^ 

0.0 

0.0 

0.0 

0.0 

50 

3.9x  10-2 

1.7x10-2 

0.1378 

0.1476 

8.4x10"* 

8.4x  10"* 

4.9x  10-'2 

4.3x  10-'2 

100 

5.98x10-° 

3.04x10-2 

6.5x10-2 

4.8x  10-2 

0.1386 

0.1430 

1.5x10-2 

1.5x10-2 

(d)  Probabilities  P(0.3l/t)  (exact)  and  P/ (0,3  Ik)  (under  IROA) 


k 

=  5 

it  = 

10 

k  = 

20 

k  = 

50 

9 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

0.999 

0.999 

0.999 

0.999 

1.0 

1.0 

I.O 

I.O 

50 

7.2x  10"* 

7.57x  10-2 

0.3467 

0.3271 

0.9975 

0.9975 

0.999 

0.999 

100 

1.14x10-2 

1.68x10-'° 

9.69x10"* 

3.72x  10"* 

0.4065 

0.3982 

0.999 

0.999 
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Table  5 


Enor  probabilities  computing  using  exact  and  IROA  models  for 
asynchronous  I^SSMA  with  RS  (32,16)  coding  and  erasures-correction  decoding 

(£,/yVo=  10  dB) 

(a)  Probabilities  P(2,0\k)  (exact)  and  P/(2,01/t)  (under  IROA) 


k  = 

=  5 

k  = 

10 

k  = 

20 

k  = 

50 

<? 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

4.6U  10"^ 

3.77x10~2 

5.94x  10"*^ 

6.01X  10"*^ 

0.0 

0.0 

0.0 

0.0 

50 

0.9999 

0.9999 

0.9863 

0.9863 

0.1565 

0.1532 

7.49x  10'*^ 

7.49x  10"‘^ 

100 

0.9999 

0.9999 

0.9999 

0.9999 

0.9797 

0.9797 

9.17x10"^ 

8.99x  10'^ 

(b)  Probabilities  P(l,l  I/:)  (exaa)  and  P/(l,l  li)  (under  IROA) 


k  = 

=  5 

k  = 

:  10 

k  = 

:  20 

k  = 

:  50 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

0.1481 

0.1565 

7.75x10"'^ 

7.75x10“'' 

0.0 

0.0 

0.0 

0.0 

50 

6.06X  10"'' 

6.06X  lO""* 

6.81x10"^ 

6.85x  10“^ 

0.2349 

0.2382 

8.65x  lO"”* 

8.65x  lO""* 

100 

2.44x10"*' 

2.44x10"" 

2.50x  10"^ 

2.50x  10“^ 

1.01x10"^ 

1.01x10"^ 

8.56x10"^ 

8.58x  10"^ 

(c)  Probabilities  P(0,2ik)  (exact)  and  P/(0,2\k)  (under  IROA) 


k  ■- 

=  5 

k  = 

10 

k  = 

20 

k  = 

50 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

0.6576 

0.6493 

0.9999 

0.9999 

1.0 

1.0 

1.0 

1.0 

50 

9.73x10"" 

3.67x  10"'^ 

8.81x10"^ 

4.75x10"^ 

0.3734 

0.3703 

0.9999 

0.9999 

100 

6.42X  lO"*'* 

5.97x  10“^^ 

7.62x10"" 

6.26x  10"'^ 

1.36x10"^ 

1.05x10-^ 

0  8195 

0.8194 
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Table  6 


Error  probabilities  computed  using  exact  and  IROA  models  for 
asynchronous  Hl/SSMA  with  RS  (32,16)  coding  and  erasures-correction  decoding 

(£^/No  =  10  dB) 

(a)  Probabilities  £(3,01/!:)  (exact)  and  Pi(3,0\k)  (under  IROA) 


k  = 

=  5 

k  = 

10 

k  = 

20 

k  = 

50 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

1.31x10"^ 

7.32x10'^ 

4.47x10”'^ 

4.66xl0~’^ 

0.9697 

0.9696 

0.0 

0.0 

50 

0.9999 

0.9999 

0.9796 

0.9795 

6.39x10-2 

5.99x10-2 

6.48X  10-*’ 

6.48X  10'” 

100 

0.9999 

0.9999 

0.9999 

0.9999 

0.9697 

0.9697 

9.06X  10-^ 

8.52X  10-^ 

(b)  Probabilities  £(2,1 1 /t)  (exact)  and  £/(2,l  l/t)  (under  IROA) 


k  = 

=  5 

k  = 

10 

k  = 

20 

k  = 

50 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

3.29x  10-2 

3.04x10-2 

5.94x1 0-‘2 

”^.04x10-'^ 

9.95x  10-2 

1.00x10-2 

0.0 

0.0 

50 

6.06x10-2 

6.06X  10-2 

6.72x  10'^ 

6.80x10-2 

9.26x  10-2 

9.32x  10-2 

7.49x  10-’2 

7.49x10-*2 

100 

2.44x10-" 

2.44x10'" 

2.50x  10"^ 

2.50x  10-^ 

9.95x10-2 

1.00x10-2 

8.27x10-2 

8.52x  10-^ 

(c)  Probabilities  P(l,Hk)  (exact)  and  £/(I,2f/!:)  (under  IROA) 


k  = 

=  5 

k  = 

10 

k  = 

20 

k  = 

<7 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

0.1152 

0.1261 

7.75x10-2 

7.75x  10'' 

0.0 

0.0 

0.0 

0.0 

50 

9.73x10'" 

3.67x10-‘2 

8.64x  10'^ 

4.72x10'^ 

0.1423 

0.1449 

8.65x  10“2 

8.65X  10'2 

100 

6.42x  10-*2 

5.97x  10-22 

7.61x10-*' 

6.26X  10-*2 

1.34x10^ 

1.03x10-^ 

7.74x  10-2 

7.77X  10-2 

(d)  Probabilities  £(0,3 1  it)  (exact)  and  £/(0.3i/t)  (under  IROA) 


k  = 

=  5 

k  = 

k  = 

k  = 

50 

<7 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

Exact 

IROA 

10 

0.5424 

0.5232 

0.9999 

0.9999 

1.0 

1.0 

1.0 

1.0 

50 

1.37x10-*2 

2.23x  10'*’ 

1.70x10-^ 

3.28x  10'2 

0.2314 

0.2254 

0.9999 

0.9999 

100 

2.78x  10-*2 

1.46x10-22 

9.58x  10"'2 

1.56x10-’2 

2.22x10-^ 

1.07x10-^ 

0.7422 

0.7417 
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Table  7 


Error  probabilities  computed  using  exact  and  Gaussian  models  for 
asynchronous  FH/SSMA  with  RS  (32,16)  coding,  Jfc=10  and  9=100. 

(a)  Probabilities  /’(2,0I^)  (exact)  and  Pg(2.01A:)  (Gaussian) 


Ei,INo  (dB) 

errors  d 

ecoding 

eiTors/erasures  decoding 

erasures  decoding 

Exact 

Gauss 

Exact 

Gauss 

Exact 

Gauss 

6 

0.6320 

0.4916 

0.9900 

0.9865 

0.9998 

0.9963 

8 

0.7782 

0.7782 

0.9999 

0.9885 

0.9999 

0.9894 

10 

0.8633 

0.7979 

0.9999 

0.9885 

0.9885 

0.9885 

12 

0.7982 

0.7982 

0.9999 

0.9885 

0.9999 

0.9885 

OO 

0.8635 

0.7982 

0.9999 

0.9885 

0.9999 

0.9885 

(b)  Probabilities  P(l,l  lik)  (exact)  and  Pc(\A\k)  (Gaussian) 


EbJNo  (dB) 

errors  decoding 

errors/erasures  decoding^ 

erasures  decoding 

Exact 

Gauss 

Exact 

Gauss 

Exact 

Gauss 

6 

0.1588 

0.2237 

4.176x10-^ 

6.728X  10-^ 

5.777X  ICr^ 

2.848x10-2 

8 

1.035x10-' 

1.034x10-’ 

3.249X  10-^ 

5.683x10-2 

3.341xl(r^ 

5.242X  10-2 

10 

6.462JC  10-2 

9.454X  10-2 

2.510x10-^ 

5.718x10-2 

2.511xl(r^ 

5.712x10-2 

12 

9.443X  10-2 

9.443X  10-2 

2.501x10'^ 

5.718x10-2 

2.501x10"^ 

5.718x10-2 

OO 

6.454x10-2 

9.443t10-2 

2.501X  10“^ 

5.718x10-2 

2.501x10^ 

5.718x10-2 

(c)  Probabilities  P  (0,2 1  k )  (exact)  and  Pq  (0,2 1  k )  (Gaussain) 


EbfNo  (dB) 

errors  decoding 

errors/erasures  decoding 

erasures  decoding  1 

Exact 

Gauss 

Exact 

Gauss 

Exact 

Gauss 

6 

"  5^018x10-^ 

5.997X  10-2 

9.901X  10-2 

4.671x10-2 

1.284x10-® 

1.175x10-^ 

8 

1.484x10-2 

1.484x10-2 

1.460x10-^ 

4.707X  10-2 

1.198x10-’° 

3.445X  10-2 

10 

7.438X  10-2 

1.296x10-2 

9.637X  10'® 

4.932x10-2 

7.660x10-" 

4.912x10-2 

12 

1.294x10-2 

1.294x10-2 

1.458x10-'° 

4.935X  10-2 

7.616x10-" 

4.934X  10-2 

OO 

7.421x10-2 

1.294x10-2 

7.616x10-" 

4.935x10-2 

7.616x10"" 

4.934X  10-2 

Table  8 


Error  probabilities  computed  using  exact  and  Gaussian  models  for 
asynchronous  FH/SSMA  with  RS  (32.16)  coding,  it=10  and  <7=100. 

(a)  Probabilities  P(3,0lJt)  (exact)  and  P^t^.OI^)  (Gaussian) 


E^INo  (dB) 

errors  decoding 

errors/eiusures  decoding 

erasures  decoding 

Exact 

Gauss 

Gauss 

Exact 

Gauss 

6 

0.5099 

0.3191 

0.9798 

0.9998 

0.9944 

8 

0.7857 

0.6876 

0.9829 

0.9999 

0.9843 

10 

0.8051 

0.7146 

0.9828 

0.9999 

0.9828 

12 

0.8054 

0.7149 

0.9828 

0.9999 

0.9828 

oo 

0.8054 

0.7149 

0.9828 

0.9999 

0.9828 

(b)  Probabilities  P(2,l  \k)  (exact)  and  Pc (2.1 1^)  (Gaussian) 


El, /No  (dB) 

errors  decoding 

errors/erasures  decoding 

erasures  decoding 

Exact 

Gauss 

Gauss 

Exact 

Gauss 

6 

0.1121 

0.1714 

6.682X  lOr^ 

5.776X  10-^ 

2.844X  10-2 

8 

6.334X  10-2 

9.059X  10-2 

5.636X  10-2 

3.340x10-^ 

3.208X  10-2 

10 

5.816x10-2 

8.333X  10-2 

5.669X  10-2 

2.511x10-* 

5.663X  10-2 

12 

5.809X  10-2 

8.324X  10-2 

5.669X  10-2 

2.501X  10-* 

5.669X  10-2 

OO 

5.809X  10-2 

8.325X  10-2 

5.669X  10-2 

2.501X  10-* 

5.669X  10-2 

(c)  Probabilities  P(1,2I/:)  (exact)  and  Pc{l,2\k)  (Gaussian) 


EbINo  (dB) 

errors  decoding 

errors/erasures  decoding 

erasures  decoding 

Exact 

Gauss 

Gauss 

Exact 

Gauss 

6 

3.678X  10-2 

5.163x10-^ 

4.637X  10-2 

1.284x10-“ 

1.171x10-* 

8 

7.67X  10-2 

1.287x10-2 

4.654X  10-2 

1.198xl0-‘° 

3.417x10-2 

10 

6.466X  10-2 

1.121x10-2 

4.873X  10-2 

7.660x10-" 

4.853x10-2 

12 

6.452X  10-^ 

1.119x10-2 

4.875X  10-2 

7.615x10-" 

4.875X  10-2 

OO 

6.452X  10-2 

1.119x10-2 

4.874X  10-2 

7.615x10-" 

4.874X  10-2 

(d)  Probabilities  P  (0,31  A:)  (exact)  and  P^  (0,3 1  A:)  (Gaussian) 


Eb/No  (dB) 

errors  decoding 

errors/erasures  decoding 

erasures  c 

lecoding 

Exact 

Gauss 

Gauss 

Exact 

Gauss 

6 

1.340x10-2 

8.183x10-^ 

3.295X  lOr^ 

6.117x10-‘2 

1.353x10-^ 

8 

1.242x10-2 

1.973x10-2 

5.319xl(r2 

7.264X  10-" 

2.741x10-2 

10 

9.718x10-^ 

1.751x10-2 

5.951x10-2 

1.169x10-‘2 

5.889X  10-2 

12 

9.687X  10-^ 

1.748x10-2 

5.959X  10-2 

8.808x1 0-'2 

5.957X  10-2 

OO 

9.687x10-^ 

1.748x10-2 

5.957X  10-2 

9.580x10-'2 

5.957X  10-2 
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Table  9 


Error  probabilities  computed  using  exact  and  Gaussian  models  for 
asynchronous  FH/SSMA  with  RS  (32,16)  coding,  EhlN^lOdB  and  q=\W. 

(a)  Probabilities  P(2,0\k)  (exact)  and  Pg (2,0 lit)  (Gaussian) 


errors  decoding 

errors/erasures  decoding 

erasures  decoding 

K 

Exact 

Gauss 

Exact 

Gauss 

Exact 

Gauss 

5 

0.9989 

0.9042 

0.9990 

0.9044 

0.9999 

0.9045 

10 

0.8635 

0.7979 

0.9152 

0.9885 

0.9999 

0.9885 

20 

7.36x  10~2 

4.0151x10-^ 

0.8921 

0.9722 

0.9797 

0.9722 

50 

2.6x  10"’° 

1.077x10"" 

7.25x  10"" 

3.891X  10"^ 

9.17x10"^ 

3.885x1 0"3 

(b)  Probabilities  P(l,l  l)k)  (exact)  and  (Gaussian) 


errors  decoding 

errors/erasures  decoding 

erasures  decoding 

K 

Exact 

Gauss 

Exact 

Gauss 

Exact 

Gauss 

5 

5.45x  10"^ 

4.5293X  10"^ 

WEgBSm 

2.44x10"" 

10 

6.45x  10~2 

9.454X  10"^ 

8.27x  10"^ 

5.718x10"^ 

■twin 

5.712x10"^ 

20 

0.1906 

1.33x10"* 

1.379x10"^ 

50 

1.50x10"^ 

4.931x10"^ 

6.247X  10"^ 

6.249X  10"^ 

(c)  Probabilities  P(0,2I/:)  (exact)  and  PG(0.2\k)  (Gaussian) 


h. 

ecoding 

erasures  decoding 

K 

Exact 

Gauss 

Exact 

Gauss 

Exact 

Gauss 

5 

6.09x  10"® 

5.170x10"^ 

8.32x  10"^ 

TnExTo^ 

6.42x10"*’ 

5.168x10"^ 

10 

7.42x10"^ 

1.296x10"^ 

3.75x10"^ 

4.932X  10"^ 

7.62x10"" 

4.912x10"^ 

20 

0.5451 

0.6366 

2.87x10"^ 

1.879X  10-^ 

1.36x10-^ 

1.867x10"^ 

50 

0.9999 

0.9999 

0.9123 

0.871 1 

0.8195 

0.8711 
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Table  10 


Error  probabilities  computed  using  exact  and  Gaussian  models  for 
asynchronous  FH/SSMA  with  RS  (32.16)  coding,  E^IN^f^lOdB  and  q=l00. 

(a)  Probabilities  P(3,0li)  (exact)  and  (3.0 lit)  (Gaussian) 


t. 

errors  decoding 

errors/erasures  decoding 

erasures  decoding 

K 

Exact 

Gauss 

Gauss 

Exact 

Gauss 

5 

0.9983 

0.8632 

0.8635 

0.9999 

0.8635 

10 

0.8054 

0.7146 

0.9828 

0.9999 

0.9828 

20 

2.10x10”^ 

7.880x10"^ 

0.9585 

0.9697 

0.9686 

50 

4.0QX 

9.0393X  10-’^ 

1.99x10'^ 

9.06X  lO"^ 

1.983x10^ 

(b)  Probabilities  P(2,l  I  A:)  (exact)  and  Pg(2.1  (Gaussian) 


l_ 

errors  c 

lecoding 

errors/erasures  decoding 

erasures  decoding 

K 

Exact 

Gauss 

Gauss 

Exact 

Gauss 

5 

5.49x  10"^ 

4.095x10"^ 

2.44x10'" 

4.083X  10'^" 

10 

8.333X  10"^ 

5.669x  10'^ 

5.664X  10'^ 

20 

9.95x  10'^ 

1.361x10'^ 

50 

8.27x10"^ 

3.687X  10'^ 

(c)  Probabilities  P(,\,2\k)  (exact)  and  Pc{l,2\k)  (Gaussian) 


errors  c 

erasures  decoding 

K 

Exaa 

Gauss 

Gauss 

Exact 

Gauss 

5 

4.336X  10'^ 

6.42x10'*’ 

mmmm 

10 

1.121x10'^ 

7.61x10'*’ 

4.853X  10'^ 

20 

0.1386 

1.34x10-^ 

50 

4.931x10-^ 

7.74x  10'^ 

(d)  Probabilities  P(0,31Jk)  (exact)  and  Pc (0.3 lit)  (Gaussian) 


t. 

errors  decoding 

erasures 

decoding 

K 

Exact 

Gauss 

Gauss 

Exact 

Gauss 

5 

1.14x10"’ 

8.346X  10-^ 

8.361x10-^ 

2.78x  10'*^ 

8.345X  10"^ 

10 

9.69X  lO"^ 

1.751x10'^ 

5.951x10'’ 

9.58x10'*’ 

5.889x10"’ 

20 

0.4065 

0.5073 

2.416x10"^ 

2.22x10"^ 

2.384x10-* 

50 

0.9990 

09999 

0.8124 

0.7422 

0.8123 
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